Preprint Number: ADP-98-033/T306 



Dimensionally Regularized Study of Nonperturbative 
Quenched QED 



Andreas W. Schreiber, Tom Sizer and Anthony G. Williams 
Special Research Centre for the Subatomic Structure of Matter 

and 

Department of Physics and Mathematical Physics, 
University of Adelaide, 5005, Australia 
E-mail: awilliam@physics. adelaide. edu. au 

We study the dimensionally regularized fermion propagator Dyson-Schwinger equa- 
tion in quenched nonperturbative QED in an arbitrary covariant gauge using the 
Curtis-Pennington vertex and perform nonperturbative renormalization numeri- 
cally. The nonperturbative fermion propagator is solved in D = 4 — 2e < 4 
dimensional Euclidean space for a large number of values of e for two values of 
the coupling, a = 0.6 and a = 1.5. Results for D = 4 are then obtained by 
extrapolation to e — > 0. We compare these results against previous studies em- 
ploying a modified ultraviolet cut-off rcgularization and find agreement to within 
the numerical precision of the present calculations. 



1 Introduction 

The Dyson-Schwinger equations (DSE) are an infinite tower of integral equa- 
tions that self-consistently relate the Green's functions of a field theory. They 
provide an alternative to lattice gauge theory El for investigating nonperturba- 
tive phenomena such as dynamical chiral symmetry breaking and dynamical 
mass generation. The advantage of such a continuum formulation is that it is at 
times possible to obtain analytical insights, there is no limit to the momentum 
range that can be studied and one is able to compare different regularization 
schemes. The disadvantage is that the tower has to be truncated to those 
equations involving Green's functions with relatively few external legs: conse- 
quently some of the Green's functions appearing in the remaining equations 
are no longer determined self-consistently by the DSE's and we must make an 
educated guess (Ansatz) regarding their form. 

As an exanTole- quantum electrodynamics (QED) in four dimensions in the 
DSE formalismoij'cl admits a dynamical mass above a certain critical coupling. 
The coupled integral equations determining the photon and fermion propaga- 
tors are completely closed once the photon-fermion proper vertex is specified. 
The quenched approximation, in which the photon propagator is bare and 
the coupling does not run with momentum, reduces the study of these coupled 
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equations to that of the fermion propagator: for computational reasons it is al- 
most universally used at present. The vertex Ansatz is constrained by_the need 
to respect the discrete symmetries and the Ward-Takahashi identities □ (WTI), 
have no artificial kinematic singularities, satisfy the requirements of perturba- 
tive multiplicative renormalizability (MR) , and agree with perturbation theory 
in the weak coupling limit. Furthermore the gauge dependence of the resulting 
fermion propagator should be consistent with the Landau-Khalatnikov trans- 
formationtj. Ultimately, observables like the critical coupling should be verified 
to be gauge-independent, once the effect of the regulator is removed. 

A number of discussions of the choice-of, P^Pf r v *? rt, fi Y ' 11 quenched QED 
can be found in the literature, e.g., Refs. D^iJ'EjIiJo'tj'ollj. Here we concen- 
trate on the Curtis-Pennington (CP) vertex, which satisfies both the WTI and 
the constraints of multiplicative renormalizability. With a bare vertex (which 
breaks both gauge invariance and MR) the critical coupling of quenched QED 
differs by approximately 50% when calculated in the Feynman and Landau 
gauge. This should be compared to a difference ofJ.css than 2% if this quantity 
is calculated with the Curtis-Pennington vertex E3E9. Similar results are ob- 
tained inta. These calculations used a form of the ultraviolet (UV) regulated 
DSE with a gauge covariance correction, where an obvious gauge covariance 
violating term arising from the cut-off was omitted. 

Clearly, the gauge dependence of the critical coupling is decreased, but not 
eliminated, through the use of a photon-fermion proper vertex which satisfies 
the WTI. The question arises whether the remaining slight gauge dependence 
in the critical coupling is primarily due to limitations of the vertex itself or 
whether it is due to the_use of a UV cut-off regulator in these calculations. 
Bashir and Pennington t§EJ have pursued the first of these alternatives and 
have obtained, within a cut-off regularized theory, further restrictions on the 
transverse part of the vertex which ensure by construction that the critical cou- 
pling indeed becomes strictly gauge independent. It is not clear that adjusting 
the vertex to remove an unwanted gauge-dependence is the most appropriate 
procedure when using a gauge-invariance (and Poincare-invariance) violating 
regularization scheme. 

We explore the second alternative by studying renormalized quenched 
QED numerically using the CP vertex and dimensional regularizationEHI where 
gauge invariance is explicitly maintained. To our knowledge this is the first 
complete nonperturbative demonstration of dimensional regularization and 
renormalization, and followsjimcfidures first developed using the DSE for a 
UV cutoff regulator in Refs.E3'li3E3. The resulting nonlinear integral equation 
for the fermion propagator is solved numerically in D = 4 — 2e < 4 Euclidean 
dimensional space. Successive calculations with decreasing e are then extrapo- 
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lated to e = and compared with the corresponding (modified and unmodified) 
UV cutoff results. 



2 Formalism 

We summarize the implementation of nonperturbative renormalization using 
dimensional regularization within the contepdjOfjiumerical DSE studies, adopt- 
ing a notation similar to that used in Refs.EO'ESIlj. The formalism is presented 
in Minkowski space and the Wick rotation into Euclidean space can then be 
performed once the equations to be solved have been written down. Note we 
use dimensional regularization in conjunction with off-shell renormalization, 
rather than using the MS or MS renormalization schemes, which can only be 
defined in a perturbative context. 

The renormalized inverse fermion propagator is defined through 

S-\fx;p) = A(v;p 2 ) jrf - B(fi;p 2 ) = Z 2 (fx, e)\j - m (e)] - e;p) 

= p'-miii) -Y,(n;p) , (1) 

where fi is the chosen renormalization scale, m(/i) is the value of the renormal- 
ized mass at p 2 — /i 2 , mo(e) is the bare mass and ^(M) e ) i s the wavefunction 
renormalization constant. Due to the WTI for the fermion-photon proper ver- 
tex, we have for the vertex renormalization constant Zi(fj,, e) = ^2(/i,e). The 
renormalized and unrcnormalized fermion self-energies are denoted as 
and E'(/^, e;p) respectively. These can be expressed in terms of Dirac and 
scalar pieces, where for example 

E'(/i,e;p) = H' d (n,e;p 2 ) p 1 +Y 1 ' s (p,e-p 2 ) , (2) 

and similarly for E(/z;p). In the interests of notational brevity we do not 
explicitly indicate the dependence on e of the renormalized quantities A(/i; j> 2 ), 
B(n;p 2 ) and £(/t;p), since for these and other renormalized quantities we 
are always interested in their e — > limit. The renormalized mass function 
M(p 2 ) = B(n;p 2 )/A(/j,;p 2 ) is renormalization point independent, which follows 
straightforwardly from multiplicative renormalizabilityE3. 
The renormalization point boundary condition 

S ~ 1 (^P)\ P 2 =I1 2 =i>- m(fi) (3) 

implies that A(p,; /i 2 ) = 1 and m{p) = M(/i 2 ) and yields the following relations 
between renormalized and unrcnormalized self-energies 

ZdAviP 2 ) = s L(/^ £ ;p 2 ) - eJuO^m 2 ) • ( 4 ) 
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Also, the wavefunction renormalization is given by 

Z 2 ( A1 ,e) = l + E' d ( M ,e;/i 2 ) (5) 

and the bare mass mo(e) is linked to the renormalized mass m(p) through 

/no (e) = [m(fjL) - £' s (/x, e; y 2 )] /Z 2 {fi, e) . (6) 

It also follows from MR that under a renormalization point transformation 
jj, — > //, m(fi') = M(p,' 2 ) and Z 2 (// , e) = A(/u'; /i 2 ) Z 2 (/it, e) as discussed inE3. 
The unrenormalized self-energy is given by the integral 

/d D k 
(2tt) d 1XS ^'' k ^^ ] k 'P) D ^{^P- k ), (7) 

where v is an arbitrary mass scale introduced in D dimensions so that the 
renormalized coupling e(/i) remains dimensionless. Since we are here working 
in the quenched approximation we have Z^{p, 1 e) = 1, eo = e(/z), and the renor- 
malized photon propagator D A " y (^; q) is equal to the bare photon propagator 

V T T J T 

with £ being the covariant gauge parameter. Finally, r l '(/x; k,p) is the renor- 
malized photon-fermion vertex for which we use the CP Ansatz, namely (q = 
k-p) 

r^ f i;k,p) = T^;k,p) + T 6 ( f i;k 2 ,p 2 ,q 2 )[r(p 2 -k 2 ) + (p + kr ( i\ , (9) 

where Tec is the usual Ball-Chiu part of the vertex which saturates the Ward- 
Takahashi identity!] 



{k+ P y 



r£ c ( M ; k,p) = ~[A(m; k 2 ) + A(p;p 2 )]r (10) 



ul [Afc k 2 ) - A{^p 2 )]L—L. - [Bfa k 2 ) - B{^p 2 )] 

and the coefficient function rg is that chosen by Curtis and Pennington, i.e., 
T 6 (n;k 2 ,p 2 ,q 2 ) = -^[A( f x;k 2 )-A( f x;p 2 )]/d(k,p), (11) 

where 

k l + pr 
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The unrcnormalizcd scalar and Dirac self-energies are extracted out of the 
DSE, Eq. (0), by taking jTt of this equation, multiplied by 1 and ft/p 2 , re- 
spectively. Note that we use the following conventionsEil for the Dirac algebra: 

9%=D => 7 m 7m - D, yy 7|t = (2 - D)^ 

Tr[l] =4 Trfy'V] = W ■ 

The integrands appearing in Eq. (Q) only depend on the magnitude of the 
internal fermion's momentum k 2 as well as the angle 9 between the fermion and 
photon momentum. Hence the D-dimensional integrals reduce to 2-dimensional 
ones, i.e., 



P P POO PIT 

d D kf(k 2 ,p 2 ,k-p) = dfl^x dkk°- x I d9 sin D " 2 6»/(fc 2 ,p 2 ,fc-p) 

(13) 

where 

is the surface area of a D-dimensional sphere. Furthermore, it is possible to 
express all the angular integrals in terms of a single hypergeometric function so 
that it is only necessary to do the momentum integral numerically, solving the 
renormalized fermion DSE by iteration. Note that the momentum integration 
extends to infinity, necessitating a change in integration variables. The infinite 
range of the integration also requires an extrapolation of A(n;p 2 ) and B(/j,;p 2 ) 
above the highest gridpoint. We check insensitivity to this extrapolation by 
comparing results obtained with a number of different extrapolation prescrip- 
tions. In addition, we use grids which extend some 20-30 orders of magnitude 
beyond what is usually used in cut-off studies. We have verified the effect of 
the extrapolation to infinity is well-controlled. 



3 Results 

We present solutions for the DSE for two values of the coupling a — e 2 /An, 
namely a — 0.6 and a — 1.5, chosen to correspond to couplings respectively 
well below and above the critical coupling found in previous UV cut-off based 
studies. Note that all results are quoted in terms of dimensionless units, i.e., all 
mass and momentum scales can be simultaneoulsy multiplied by any desired 
mass scale. 

Fig. |l| show a family of solutions with the regulator parameter e decreased 
from 0.08 to 0.03 for the two values of the coupling. We see that the mass 
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0t=0.6 solutions for different e 



a=1.5 solutions for different e 




Figure 1: The finite renormalization A(p 2 ) and mass function M(p 2 ) for various choices of 
the regulator parameter e. Figure 1(a) has coupling a = 0.6 and 1(b) a = 1.5, with gauge 
parameter £ = 0.25, renormalization point fi 2 = 10 s , renormalized mass m(fj,) = 400 and 
scale v = 1. In the low p 2 region the smallest e has the largest value of M(p 2 ). 



function increases in strength in the infrared and tails off faster in the ultravi- 
olet as e is reduced or a is increased. Furthermore, it is important to note the 
strong dependence on e, even though this parameter is already rather small. 
As one would expect, the ultraviolet is most sensitive to this regulator, how- 
ever even in the infrared there is considerable dependence due to the intrinsic 
coupling between these regions by the renormalization procedure. This strong 
dependence on e should be contrasted with the situation in cut-off based stud- 
ies where it was observed that already at rather modest cut-offs (A 2 ks 10 10 ) 
the renormalized functions A and M had reached their asymptotic limits. 

Currently it is not possible to decrease e significantly below the values 
shown in Fig. [j] due to numerical limitations. In order to extract the values 
of A and M in four dimensions we extrapolate to e = 0, which introduces 
an added uncertainty in the final result. However it is possible to estimate 
this uncertainty by using the fact that in the limit e — > the renormalized 
quantities should become independent of the arbitrary scale introduced to 
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keep the coupling a dimensionless in D dimensions. 




Figure 2: (a) shows the finite renormalization A(p 2 ) and mass function M(p 2 ) evaluated at 
p 2 = 1 for various values of the regulator parameter t and extrapolated to e = by fitting 
a polynomial cubic in e. Shown are five scales v = 1, 10, 100, 1000 and 10000. All other 
parameters are those of Fig. 1. The different scales coincide at e = to an accuracy of 
approximately 0.2%. (b) shows the finite renormalization A(p 2 ) and mass function M(p 2 ) 
extrapolated to e = at each momentum point for the three scales v = 100, 1000 and 10000. 
Also shown are results obtained using a modified and unmodified UV cutoff. The modified 
UV cutoff and extrapolated results are indistinguishable in these plots, and differ from those 
of the unmodified UV cutoff in the IR. 

In Fig. H^, we show A(p 2 ) and M(p 2 ) evaluated with a = e^/Air = 1.5 in 
the infrared (at p 2 = 1) as a function of e for a range of values of v. The 
results at e = are extracted from cubic polynomial fits in e. The agreement 
between the different curves at e = is excellent, being of the order of 0.2%. 
In Fig. we show the results for v = 100, 1000 and 10000 extrapolated to 
e = as a function of momentum for a = 1.5. Again the agreement is very 
good for a wide range of p 2 and v. Only in the ultraviolet region (above say 
p 2 = 10 12 ) are differences between the curves discernible. The discrepancies 
for v — 1 and 10 cases were somewhat greater in the UV, presumably reflecting 
greater non-linearity in the fits for these cases that is manifest at p 2 = in 
Fig. ||a. This is supported by the observation that these lower-scale cases were 
more sensitive to changing the order of the polynomial the higher-scale cases. 
Also shown in Fig. Hb are the corresponding UV cutoff results, both with and 
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without a modification which used a gauge-covariarice fix to remove an obvious 
part of the gauge dependence induced by the cutoffEa. The modified UV cutoff 
curve is indistinguishable from the scaled dimensionally regularized ones, while 
the unmodified UV cutoff curve clearly deviates from the others in the IR. It 
should alsp_be noted that the oscillatory behaviour in the mass function first 
noticed in EH is reproduced in this work. 

4 Conclusions and Outlook 

We have reported the first detailed study of the numerical renormalization 
of the fermion Dyson-Schwinger equation of QED through the use of a di- 
mensional regulator rather than a gauge invariance- violating UV cut-off. The 
initial results are encouraging. Firstly, we have explicitly demonstrated that 
the approach works and is independent of the intermediate dimensional reg- 
ularization scale as expected. Secondly, the previously obtained modified UV 
cut-off results are reproduced with considerable precision. 

A significant practical difference between the dimensionally and UV cut- 
off regularized approaches is that in the former it is currently necessary to 
perform an explicit extrapolation to e = whereas in the latter for a sufficiently 
large choice of UV cut-off the results became independent of the cut-off. We 
are presently investigating whether it is numerically possible to extend the 
studies to smaller values of e in order to improve the precision of the e — * 
extrapolation. In the meantime, the extrapolation to e = appears to be well 
under control, at least for values of the fermion momentum away from the 
ultraviolet region. 

Having demonstrated the numerical procedure of renormalization using 
dimensional regularization, we plan to study chiral symmetry breaking and 
in particular hope to extract the critical coupling as a function of the gauge 
parameter, and to shed further light on the nature of the chiral limit discussed 
in til Results of this work will be presented elsewhere. The eventual aim is to 
extend this treatment to the case of unquenched QED. 
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